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Abstract 

We apply the Green function formalism for t — t production and de- 
cay near threshold in a study of the effects due to the momentum 
dependent width for such a system. We point out that these effects 
are likely to be much smaller than expected from the reduction of 
the available phase space. The Lippmann-Schwinger equation for the 
QCD chromostatic potential is solved numerically for S partial wave. 
We compare the results on the total cross section, top quark intrinsic 
momentum distributions and on the energy spectra of W bosons from 
top quark decays with those obtained for the constant width. 
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1 Introduction 



There are accumulating indications that the long awaited top quark will be 
discovered soon. The value of its mass m t preferred by the indirect LEP 
measurements Q implies that the t quark is a short-lived particle, and its 
width F t is of the order of several hundreds MeV. As a consequence the cross 
section for ti pair production near energy threshold has a rather simple and 
smooth shape. In particular, it is likely that in e + e~ annihilation only the 15" 
peak survives as a remnant of toponium resonances [§, § . Nevertheless, as 
pointed out by Fadin and Khoze || the excitation curve cx(e + e~ — > ti) allows 
a precise determination of m t as well as of other physical quantities such as T t 
and the strong coupling constant a s . These results, derived analytically for 
the Coulomb chromostatic potential, have been confirmed by Strassler and 
Peskin || who studied numerically a more realistic QCD potential. The idea 
[|], [l| to use Green function instead of summing over overlapping resonances 
has been also applied in calculations of differential cross sections: energy 
spectra of W's from decays || , and intrinsic momentum distributions of top 
quarks in ti systems || [7| . 

Pandora's box was opened when the authors of @ found that the effects 
of the momentum dependent width may show up in the annihilation cross 
section cr(e + e~ — > ti) . Despite the fact that the calculation is plagued 
by unpleasant and unphysical gauge dependence, one has to consider their 
conclusions as an indication that the accuracy of the calculations ^ || 
assuming a momentum independent decay rate should be reconsidered. One 
cannot use the results of 0] for such an estimation, because some potentially 
relevant effects are not taken into account there. 

In the present paper we extend the calculation of by including a mo- 
mentum dependent width. Our work is based on an observation M that 
relatively simple semiclassical arguments reproduce well the results of a com- 
plete calculation of the lifetimes and the energy spectra for decays of muons 
bound in light and medium nuclei. Since the strength of interaction is com- 
parable for ti and /i~-nucleus for the nucleus charge Z ~ 10-20, one may 
expect that analogous considerations closely approximate the results of the 
complete QCD calculations for top quark pair production near threshold. 

One of the most important future applications of this calculation is the 
determination of m t and a s from the combined measurements of the total and 
the differential cross sections in e + e~ — > ti . We study the theoretical uncer- 
tainties related to model assumptions on the momentum dependent width. 
These result from 0(a^) corrections to the process e + e~ — > ti — > bW + bW~ 
which have been not calculated yet. Some corrections due to reduction of the 
available phase space || have been considered in |7|]. However, as we already 
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pointed out in our previous paper || one should expect large cancellations 
between the corrections due to the phase space reduction and the Coulomb 
enhancement of the b and b wave functions. In order to determine theoretical 
uncertainties related to the momentum dependent width we evaluate predic- 
tions of a model (Model I) which overestimates the effects on a(e + e~ — > ti) 
and the differential cross sections. Nevertheless, the resulting theoretical un- 
certainties in determination of m t and a s are quite small. 
We also work out the predictions of a more realistic model (Model II) for the 
momentum dependent width and find out that for this model the results are 
close to those obtained assuming constant width. 

Our paper is organized as follows: In sect. 2 we describe the physics of 
bound ti decays and formulate our models of the momentum dependent 
width. In sect. 3 , which we add for the sake of completeness, the Green 
function method and our QCD potential are briefly described. In sect. 4 the 
results of the numerical calculations are presented, and in sect. 5 our conclu- 
sions are given. 

2 Models for the width of t — t system 

In our earlier paper |J we followed |5| and confined the discussion to 
the non-relativistic approximation and the constant decay rate. The basic 
reason was that numerous effects related to the intrinsic momentum in the 
t—i system were known to be small and they were contributing with opposite 
signs to the width. Let us list now the main physical effects which result in 
the momentum dependent width. The width of the t — i system depends on 
the intrinsic momentum, of say t quark, because both the matrix element and 
the phase space available for the decay products depend on it. The phase 
space effect tends to reduce the decay rate of bound top quarks relative to 
free ones ||, and the effect is enhanced, because for short-lived particles the 
distribution of intrinsic momentum is broad. However, for the same reason 
the decays take place at short relative distances, where the wave functions 
of b and b quarks originating from the decays are distorted (enhanced) by 
Coulomb attraction. Therefore, when calculating the amplitude of t — > bW 
transition, one should use Coulomb wave functions rather than plane waves 
for b quarks. This effect clearly increases the rate. A third factor is the time 
dilatation: a top quark moving with velocity v lives longer in the center-of- 
mass laboratory frame]]. Let us consider now QCD corrections to the decay 

3 In an analogous calculation for fi~ bound in a nucleus of charge Z one obtains [^| 

r = r /ree [1 - 5(Zaf] [1 + 5(Za) 2 ] [l - (Zaf/2] 

where the first correction factor comes from the phase space suppression, the second from 
the Coulomb enhancement, and the third one from time dilatation. 
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t — > bW . The calculations of the corrections to the total rate |H| , the energy 



spectrum of the charged lepton ]TT] and the energy of W |12[ do not include 
two effects which may be non-negligible for t — i production near threshold: 
a) the QCD Chudakov effect, i.e. the suppression of soft gluon emission from 
a small and colour singlet ti source, and b) the additional suppression of 
the transition from a colour singlet into a colour octet ti system due to the 
related change of its potential energy^. In the following we ignore all these 
effects, as well as analogous effects in transverse gluon QCD correction to the 
production vertex. Since 0(a s ) QCD corrections to the decay JT0|, [11], [12], [13 
are ~10%, we hope that the accuracy of our approximation is of the order 
of 1%, i.e. that it is comparable to the size of a 2 effects which are also 
neglected. In particular we neglect all relativistic and spin corrections but 
the time dilatation in the momentum dependent width. QCD corrections to 
the decay rate T t _f are included as an overall factor / modifying the width 
T t for single top decay: 



T(p,E) 
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C(p,E)T t 
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f(y) = n 2 + 2Li 2 (y) ■ 
(5 + Ay) ln(l 



2Li 2 (l - y) + [Ay(l ■ 
■y)-{l-y){5 + 9y 



y-2y 2 )\ny + 2(l-y) 2 - 
■6y 2 ))/[2(l-yf(l+2y)) (3) 



The correction factor C(p, E) depends on the intrinsic momentum p and 
the energy E = y/s — 2m t as explained at the beginning of this section. 
Some effects such as the phase space reduction or the time dilatation can 
be easily implemented. In fact these effects have been already included in 
the calculations performed in 0. The Coulomb enhancement cannot be 
easily taken into account. In principle one has to replace the plane wave 
functions for b quarks by relativistic Coulomb functions when evaluating the 
amplitude for t bW transition. One may hope, however, that the following 
observation, valid for muons bound in nuclei 0, holds also for chromostatic 
attraction in t — i systems: the phase space suppression and the Coulomb 
enhancement nearly cancel each other. For light nuclei the result is well 
described by the time dilatation suppression alone. For heavier nuclei there 
is a further reduction in the decay rate because some of the decay electrons 
will not have sufficient kinetic energy to escape to infinity from the region 
of the electrostatic potential of the nucleus. We think that for t — i it is 



4 For a recent discussion of the effect of T t on soft gluon radiation see jl4| 
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reasonable to replace this 'escape into infinity' condition by a condition on 
the effective mass \x h i of the b — b system resulting from the decays of t and 
t. 

At first we discuss a Model I which overestimates the effects of the momentum 
dependent width. Comparing the suppression factors corresponding to Model 
I with the ones discussed in J7|, one observes that they are not very different, 
see Fig. |l]b in the present paper and Fig. 6 in ]7|. Thus, the results derived 
from Model I should resemble the ones obtained in 0. We assume that for 
a t — t system the t quark with the intrinsic three-momentum p decays in 
its rest frame as a free particle of mass m t . W + originates from this decay 
with energy fixed by two-body kinematics. Then its four-momentum in the 
t — i rest frame is obtained from the Lorentz transformation. The same 
procedure is applied to W~ from the decay of i. Ws from the decays are 
colour singlets, so they escape without any final state interaction. On the 
other hand the b and b interact and their momenta are changed by final state 
interactions. The four-momentum of the b — b system is fixed by energy- 
momentum conservation. It is clear that for such a model the energy spectra 
of Ws are broadened due to the Fermi motion. The effective mass can be 
smaller than for decays of unbound t and t. For some configurations of the 
momenta of W's the resulting fSr can even be negative. These configurations 
must be rejected. We impose a stronger condition on p, b i and require that \x h i 
must be larger than /i = 2rrib + A , where A= 2 GeV. It is plausible that if 
this condition is fulfilled the b — b system decays into hadrons with probability 
one. For fi b i < p, we assume complete suppression. Thus, our Model I for 
the momentum dependent width is defined by the following formula: 

dr(p, E) = T t ^l-p/(p + m 2) 6 ( m - fio) dL 2 ® 2 (4) 

where the volume of the phase space Z/2®2 for ibW + ){bW~) is normalized to 
one. Thus, if the suppression factors were absent the width T t _i of the t — t 
system would be equal to the sum of the widths of free t and t. The first 
suppression factor in (|) is due to the time dilatation, and the step 9 function 
describes the additional suppression. Depending on the total energy and the 
intrinsic momentum of the decaying top quarks, the above-mentioned effects 
result in a decrease of the width T t _i. In Figs. |l]a-b the correction factors 
C(p, E) are plotted for the top masses 120 and 150 GeV, respectively, and 
for three energies near the threshold. The phase space reduction is larger 
for energies below threshold and leads together with the relativistic time 
dilatation factor to a strong suppression of decays for large values of the 
intrinsic momentum. Both effects also depend on the top mass: the larger 
m t , the smaller the reduction of T t _f. The energy dependence is also weaker 
for smaller values of A . However, the latter effect is rather small. 
Let us describe now Model II which we consider to be more realistic. Instead 
of @) the momentum dependent width is given by the formula: 
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Figure 1: The suppression factors C(p, E) for the decays of t — i system in 
comparison to the single quark decays of t and i: a) mt = 120 GeV and b) 
m t = 150 GeV. Model I: the energies correspond to IS peak — dashed, E = 
- dotted, and E = 2 GeV — dashed-dotted lines. Model II: energies of IS 
peak — solid lines. 



dr(p, E) = Y t yjl - + (j?) 6 ( m - /i ) dL 2 ® 2 (5) 

where 

H 2 = PaP a 

and p a is the four-momentum of the off-shell top quark. In its rest frame the 
energy of W originating from the decay is given by the two-body kinematics: 

fx 2 + m 2 w - m 2 b 
Ew = 2^ (6) 

The correction factor C(p, E) evaluated from Model II depends weakly on 
energy and is close to 1, see solid lines in Figs. [l|a-b. Therefore , the corre- 
sponding results are much closer to those derived assuming constant width. 
Two more remarks are in order here. 

The models which we consider are inspired by the observation |1| that for 
short-lived top quarks the final states hadrons originate from two-jet b — b 
configurations. The lower the value of top mass and the longer its lifetime, 
the less realistic become our assumptions. For small T t the decays take place 
in a well-defined sequence, and the treatment of the recoil effects given in [Q 
is more appropriate. 



The width for off-shell top quarks is a gauge dependent object. The re- 
sults for physical quantities do not depend on the choice of gauge because the 
gauge dependent contributions from propagators, vertex corrections, wave 
function renormalization and non-resonant graphs cancel order by order in 
perturbation theory. The problem of gauge dependence in Z° and W pair 
production has been recently discussed in the literature [TEfl . The case of 



tt production near threshold is even more complicated because of multiple 
gluon exchange and resonance formation due to Coulomb-like strong attrac- 
tion. In our opinion the only way to avoid gauge dependence in final results is 
to solve the bound state problem including the constant on-shell width, and 
to treat the difference between on-shell and off-shell widths in the frame- 
work of perturbation theory. However, such a complete approach is beyond 
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the scope of the present paper. Our main goal is to estimate the importance 
of the effects related to the momentum dependent width. Therefore, we 
follow JFj and include the momentum dependent width but ignore all other 
contributions. We argue that due to the cancellations between phase space 
suppression and Coulomb enhancement the difference between the on-shell 
and off-shell widths is strongly reduced. Thus, the related differences in the 
cross sections are small, and their gauge dependence may be unimportant 
from the practical point of view. 



3 Green function for t — t system 

In this section we briefly describe the Green function method for e + e~ — > ti 
annihilation and the numerical solution of the Lippmann-Schwinger equation 
near the energy threshold; see || for details. In our discussion we neglect Z° 
contribution and transverse gluon correction to the production vertex. The 
differential cross section for top quark pair production in electron positron 
annihilation reads f\ 



^-(p,E) = 3 -^}T( P ,E)\g(p,E)f . (7) 
a A p 7r s mf 

The Green function Q(p, E) is the solution of the Lippmann-Schwinger equa- 
tion 

g(p, e) = g (p, e) + g (p, e) J S^v(p- q)g{q, e) (8) 

where V(p) is the potential in momentum space. The free Hamiltonian that 
is used to define the Green function (?o includes the momentum dependent 
width: 

So(p,E)= * (9) 

E - L + lF (p, E ) 



mi 



In our numerical calculations we use the two-loop QCD potential |L7j for 
n/ = 5 quark flavours, and join it smoothly to a Richardson-like phenomeno- 
logical potential |I§ for intermediate and small momenta. This potential 



gives a very good description of J /if) and T families. 
We define the potential V(p), where p = \p\ as follows: 



VW P ) = -^^# + Vo*) (10) 

6 p z 



where 



(Xpert 

(p) if p > Pl , 

".//(/') = { <*r(p) + (P - p 2 ) ^ik^nM i { p 1>p> p 2 , (n) 

a R (p) if p < p 2 • 



5 In (7) of H a factor (2tt) 3 is omitted. The results, plots and conclusions presented 
there are not affected by this omission. 



6 



Figure 2: cn e //(g) for different values of a s (Mz): solid: 0.12, dashed: 0.11, 
dashed- dotted: 0.13, dotted: 0.10 and 0.14- 



Figure 3: QCD-Potential in the position space Vjkt{^) for different values 
of a s (Mz) : solid: 0.12, dashed: 0.11, dashed- dotted: 0.13, dotted: 0.10 and 
0.14. 



with 
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a 
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The Richardson-like part is given by 



<xr(p) 



Ait 
~9 



1 



Cltt 



(14) 



ln(l+pVA|) P 2 (P 2 + ?L). 
Throughout our calculations we take pi = 5 GeV, p 2 = 2 GeV, A# = 400 
MeV and q cu t = 50 MeV . The constant Vo is fixed by the requirement that 
after Fourier transform the potential in the position space 

V JKT (r = lGeV" 1 ) = -l/4GeV. 

In Fig. ^|we plot a e ff for a s (Mz) varying between 0.10 and 0.14 whereas in 
Fig. |3| our potential VjKTij") in the position space is shown. 
Near the energy threshold one can neglect all but S partial waves and nu- 
merically solve the corresponding one-dimensional integral equation. The 
spherically symmetric solution fulfills the unitarity condition 

d 3 p 



T(p,E) \G(p,E)\ 2 = -ImG(f = 0,f = 0,E) 



(15) 



which for the constant decay rate reduces to the formula for the total cross 
section derived in H. 
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Figure 4: Comparison of the differential cross sections ^ evaluated for the 
momentum dependent T(p, E) (Model I: dotted, Model II: solid) and constant 
T t widths (dashed lines) for m t = 120 GeV; and a)E = —2.3, b)E = and 
c)E = 2 GeV, respectively. 



Figure 5: Comparison of the annihilation cross sections a(e + e~ — > tt ) (in 
units of R) evaluated for the momentum dependent T(p,E) (Model I: dotted, 
Model II: solid) and constant T t widths (dashed lines). 

4 Results 

As desribed in the last section, the momentum dependent width T(p, E) 
enters our numerical calculation of the Green function Q(p,E). The differ- 
ential cross sections ^ derived from (|7j) are different from the analogous 
cross sections obtained assuming constant width. In the latter case is 

° It dp 

proportional to \pG(p, E)\ 2 , and we show these distributions as dashed lines 
in Figs, ^a-c for mt = 120 GeV and T t = 0.3 GeV . For a comparison we 
plot also analogously normalized cross sections r ^ £ - > \pG(p, E)\ 2 for the mo- 
mentum dependent width (dotted lines for Model I and solid for Model II). 
The main effect of the varying width T(p, E) in Model I is a change in the 
normalization of the distribution 15% for the IS" peak) accompanied by 
the expected suppression of the large momentum tail. Thus the Green func- 
tions are slightly shifted to smaller momenta. However, it is remarkable that 
the position of the maximum is not much affected. The energy dependence 
of the annihilation cross section in the threshold region can be obtained by 
integration of the differential cross section over the intrinsic momentum: 

1 firv^ POO 

o-(e + e~ ^tt)\ E = — / dp P 2 T( P , E)\G(p, E)\ 2 (16) 
Ssrrif Jo 

In Figs. |5]a-b we show cr(e + e~ — > ti) as functions of energy for two values of 
mt '■ 120 and 150 GeV. Once again we observe a stronger effect of the mo- 
mentum dependent width for smaller m t : while for a top mass of 120 GeV 
an effective decrease of the width in Model I leads to the increase of the cross 
section at 15 resonance by more than 15%, the total cross section is much 
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Figure 6: Normalized energy spectra ofW bosons originating from the decays 
of t — t systems at IS peaks (narrower) and for E = 2 GeV (broader). The 
dashed lines correspond to the momentum dependent width evaluated in Model 

I and solid ones to the constant widths. The W spectra obtained from Model 

II are very close to the solid curves (constant width case). 



less affected for top masses as high as 150 GeV and larger. Above the 15* 
peaks, and in particular above the threshold, the effects of the momentum 
dependent width are small, typically of the order of few percent. We observe 
that also for the total cross section the position of the IS peak is not shifted, 
despite the fact that in Model I the normalization is significantly changed. 
It can be shown [BJ| that the position of the peaks in the total a(e + e~ — > ti) 



and differential ^ cross sections are essential for a precise simultaneous de- 
termination of m t and a s from the threshold region. Thus we conclude that 
the precision of such an analysis is not affected by theoretical uncertainties 
related to the momentum dependent width. 



As expected, the results obtained from Model II are close (within a few 
percent) to the results corresponding to the constant width, c.f. solid and 
dashed lines in Fig. (| and |5|. 

The top momentum distributions themselves cannot be measured exper- 
imentally. They can be reconstructed kinematically from the four-momenta 
of the final state particles, but, as should be clear from our discussion, those 
are sensitive to model assumptions. In particular the four-momenta of b 
and b quarks are distorted by final state interactions. The final state interac- 
tions, however, do not change the four-momenta of W's which are colourless. 
Thus their distributions contain important information on the dynamics of 
production and decay of top quarks. In Figs. |6|a-b we compare the normal- 
ized energy spectra of W's calculated for the momentum dependent T(p, E) 
in Model I and constant T t widths. The calculations have been performed 
for the non-zero width of W bosons, Y w = 2.2 GeV. Since the distributions 
of the intrinsic momentum are narrower for T(p,E), c.f. Figs. ^ and [5], the 
resulting energy spectra of W's are also narrower than the corresponding 
distributions for the constant width. We show these distributions as dashed 
and solid lines for the varying and constant width, respectively. The sharper 
spectra are obtained for the energies corresponding to 15 peaks and the 
broader ones for E = 2 GeV. The spectra of W's obtained from Model II are 
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close to those for the constant width. 

It it self-evident that for narrow resonances the effects related to changes 
of the width become enhanced. Thus the observed mass dependence is under- 
standable. It is clear also that the initial state radiation, not included in our 
calculations, dilutes sharp resonances and significantly reduces effects due 
to changes of the width. Therefore, in our opinion, for m t around 150 GeV 
and larger theoretical uncertainties introduced by model assumptions of the 
present paper, or those introduced in J/J , may be neglected. In particular the 
position of IS peak is rather insensitive to model assumptions on the width. 
The normalization of the cross section at the peak is sensitive to these as- 
sumptions if the top mass is in the lower range of the allowed region. Thus a 
rigorous (not model) calculation including the effects discussed in this paper 
may prove to be indispensable if the top mass is around 120 GeV and one 
aims at a high precision study based on the absolute normalization of the 
cross section at IS" peak. 

However, we would like to stress once again that considering only phase 
space suppression one overestimates effects of the momentum dependent 
width. It is likely that rigorous calculations will give results similar to those 
of Model II, which are quite close to the predictions obtained using the con- 
stant width. 

5 Summary 

In this paper we have studied the influence of the energy and momentum 
dependent width T(p, E) on t — t production near the energy threshold. The 
position of the IS* peak in the total cross section as well as the maximum of 
^ are rather insensitive to model assumptions related to T(p, E). The total 
cross section <r(e + e~ — > ti) around the IS peak reflects the effective change 
of the width introduced by T(p, E). The effects of the momentum dependent 
width are likely to be much smaller than expected from the arguments based 
on the reduction of the available phase space. Further dilution of the IS peak 
due to initial state radiation results in additional reduction of the differences 
between the results for the momentum dependent and constant widths. The 
energy spectra of W bosons are slightly narrower, because the momentum 
dependence in T(p, E) results in suppression of the large p components of 
t — i wave functions. 
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